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Abstract We address the issue of nuclear attenuation in 
nucleon and pion knockout reactions. A selection of re- 
sults from a model based on a relativistic multiple-scattering 
approximation is presented. We show transparency calcu- 
lations for pion electroproduction on several nuclei, where 
data are in very good agreement with calculations includ- 
ing color transparency. Secondly, we discuss the density 
dependence of reactions involving one or double proton 
knockout. The latter reaction succeeds in probing the high 
density regions in the deep interior of the nucleus. 
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1 Introduction 

It is an open issue how and at what energy scale hadrons 
emerge from quarks and gluons, the fundamental con- 
stituents of quantum chromodynamics. The behavior of 
nucleons in a nucleus is another topic of great interest 
and this can be studied by exploring the limits of the shell- 
model description of nuclei. This can be achieved, for in- 
stance, by looking at the medium modifications of hadron 
properties, the density dependence of the nucleon-nucleon 
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interactions, or the properties of the nucleon-nucleon cor- 
relations at very short internucleon distances. 



These fascinating issues can be experimentally explored 
at facilities that can probe the nucleonic and sub-nucleonic 
length scales, such as Jefferson lab, J-PARC, and FAIR. 
Various aspects of nuclear structure can be studied in re- 
actions such as A{e,e'p), A{e,e' ti), A(p,2p), A{j,2p), etc., 
whereby a fast leptonic or hadronic probe ejects one or 
more hadrons from the target nucleus. To interpret the 
data from these experiments, one needs models that can 
quantify the effect of nuclear attenuation on the imping- 
ing and ejected nucleons and pions. These model calcu- 
lations can for example be used to identify QCD -mediated 
deviations from traditional nuclear physics predictions and 
to map the density regions of nuclei contributing to the 
cross section of some specific reaction. 



Here, we present a selection of results obtained in a 
model based on a relativistic extension of Glauber multiple- 
scattering theory for the description of nuclear attenua- 
tion on the impinging hadrons (initial-state interactions, 
ISI) or ejected hadrons (final-state interactions, FSI). This 
relativistic multiple-scatttering Glauber approximation (RMSGA) 
model provides a comprehensive theoretical framework 
that can be used for a variety of leptonic and hadronic 
probes and one or more outgoing nucleons and/or pions. 
The RMSGA model has no free parameters. It is used to 
calculate several observables like cross sections and nu- 
clear transparencies. The general features of the model 
are introduced in Sec. [2] and some results are shown in 
Sec.|3j 



2 Relativistic multiple-scattering Glauber 
approximation 

Glauber multiple-scattering theory can be applied when 
the wavelength of the scattering particle is small in com- 
parison with the typical interaction range between the 
scatterer and the spectator particles. Examples of such 
kinematic conditions will be discussed in Sec. pj] with nu- 
cleon and pion energies in the range of a few GeV Origi- 
nating from optics, Glauber theory uses the eikonal ap- 
proximation for high energy particles scattering under 
small angles. A relativistic extension of this eikonal ap- 
proximations was developed and has been applied to A(e,e'p) 
reactions [1,2], A(p,2p) processes [3]|4]], pion photo- and 
electroproduction reactions [5,6] and neutrino-induced 
nucleon knockout 0. Relativity is accommodated both 
in the kinematics and dynamics, and all hadrons involved 
in the reaction process are described by relativistic wave 
functions. The impulse approximation is used for the in- 
teraction of the beam particle with the target. 

As an example we mention the Glauber phase in the 
RMSGA model entering the amplitude for a reaction with 
an impinging proton which has i hadrons in the final state 
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Here, a stands for the quantum numbers of the nucleon 
that interacts with the initial probe, r = (b,z) is the coor- 
dinate of the hard interaction, with z along the momen- 
tum of the impinging hadron. The coordinate r, = (b;,z,) 
is defined such that the z ; axis is along the direction of 
the momentum of the ejected hadron i. The r' { = (b^,z[) 
denotes the coordinate r' in the reference system defined 
byr,-. 

The product a occ runs over all residual nucleons, with 
4> a (V) the single particle bound-state wave functions 
obtained from the Serot-Walecka model []8]. The Heavi- 
side function ensures that only nucleons in the backward 
path of the impinging hadron and nucleons in the for- 
ward path of the ejected hadrons can make contributions 
to their phase. Information on the ISI and FSI is contained 
in the profile function, and the product i runs over all par- 
ticles subject to nuclear attenuations. The profile function 
r iN for nucleon-nucleon and pion-nucleon scattering de- 
pends on three energy-dependent parameters and has a 
Gaussian form: 
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Fig. 1 : The effect of the lower components in the wave 
functions for the scattering centers on the computed 
Glauber phase ^ at p = 1.5 GeV for proton emission from 
12 C. The figure shows the ratio of the norm of the Glauber 
phase as computed with the full relativistic single-particle 
wave function to the one which only retains the upper 
Dirac components. The ^(b,z) is computed for a proton 
escaping along the z axis. The coordinates (r, 6) indicate 
the position of the initial absorption process which trig- 
gers the reaction. 



(2) 



The three parameters cr™,e iN ,l3 iN , are fitted to nucleon- 
nucleon and pion-nucleon scattering data H9l fToUTTH . 

The scalar Glauber phase ^(b,z) is a complex quan- 
tity which encodes the combined effect of attenuation for 
a specific exclusive reaction that is initiated at the posi- 
tion r(b,z). The deviation of the norm | ^(b,z) | from 1 is 
a measure of the magnitude of nuclear attenuation. 

In our relativistic formulation of Glauber multiple-scattering 
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in Eq. (1) 



(withi = n or N'~). 



theory the single-particle densities 
contain an upper and lower component. It is conjectured 
that the effect of relativity on the magnitude of the nu- 
clear attenuation can be estimated by the relative contri- 
bution from the lower components. Fig.[T]shows the effect 
of the lower components of the single particle densities on 
the Glauber phase of an outgoing particle with an energy 
representative for the few GeV range. This shows that the 
effect of relativity on the computed magnitude of nuclear 
attenuation is rather small, typically a few percent. 

It can be expected that the nucleon-nucleon and pion- 
nucleon profile functions entering Eq. ([I]) will be sub- 
ject to medium modifications II12II . Mechanisms like Pauli 
blocking often lead to an effective reduction of the nucleon- 
nucleon cross sections in the medium. At higher ener- 
gies the effect of Pauli blocking is small. Another impor- 



tant source of medium effects are short-range correla- 
tions. Short-range correlations (SRC) can be included in 
the treatment of the Glauber eikonal phase by using the 
information that a nucleon is present at the coordinate of 
the hard interaction. 

This can be achieved in the following way. First, the 
squared single-particle wave functions in Eq. ([!]) can be 
approximated by the one-body density of the target nu- 
cleus p A ' (r) defined as 
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Here, tf^' s ' is the ground-state wave function of the target 
nucleus, obtained by antisymmetrizing the product of the 
single-particle wave functions <fi a . This substitution has 
minor effects on the results obtained in the RMSGA model 
[U . In a next step , the averaged density p A J (r) can be 
substituted with the ratio of the two-body density p A J 
(normalized as J dr 1 J dr 2 p A (r 1 ,r 2 )=A(A— 1)) and the 
one-body density: 
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where r is the coordinate of the hard interaction. In the 
case of an uncorrelated two-body density: 
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and Eq. ([3]) becomes trivial. One can include SRC in the 
two-body density by adopting the following functional 
form lEGOl : 
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where g{r 17 ) is the so-called Jastrow correlation func- 
tion and j{r) a function that imposes the normalization 
of the two-body density obtained as the solution of an 
integral equation. With the above expression for the two- 
body density , Eq. ([3]) becomes 
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P A J (r 2 ) - T(r 2 )p A iJ (r 2 )r(r)g(|r 2 -r|) = p A tt (r 2 ,r). (6) 

With above derivation it is clear that the calculation of the 
FSI effects can be corrected for SRC by replacing 4> a (V) 

with pf f (r',r)M in Eq. ^ 

Colour transparency is a QCD mediated phenomenon 
that predicts reduced hadron-nucleon interactions that 



become more pronounced as the initial hadron produc- 
tion process occurs at smaller and smaller scales. The ef- 
fect of colour transparency can be included in the RMSGA 
calculations by making use of the quantum diffusion model 
of Ref. [14] and replacing the total cross section param- 
eter in Eq. ^ with an effective one that evolves from 
a reduced to its normal value along a certain formation 
length l h : 
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Here, n is the number of elementary fields (2 for the 
pion, 3 for the nucleon), k t = 0.350 GeV/c is the average 
transverse momentum of a quark inside a hadron, and 
$€ is the hard-scale parameter (or virtuality) that gov- 
erns the CT effect. For the formation length l h k 2p/AM 2 , 
we adopt the values AM 2 = 1 GeV 2 for the proton and 
AM 2 = 0.7 GeV 2 for the pion. 



3 Results 

Transparencies First we show some results for transparency 
calculations in the RMSGA model. The transparency ob- 
servable is used in the search for the emergence of par- 
tonic degrees of freedom in nuclear reactions and is de- 
fined as the ratio of the cross section on a nucleus to A 
times the cross section on a free proton. Figure [2] presents 
the results from our transparency calculations for the pion 
electroproduction (A(e,e / 7i + )) reaction on four different 
nuclei, together with the experimental data [16] and re- 
sults from the semiclassical model of Ref. [15] . The RMSGA 
and semiclassical transparencies display a modest increase 
over the Q 2 range. This behaviour finds a simple explana- 
tion in the p n dependence of the o- t0 + . The results con- 



(5) tained in Fig. [2] cover a range in pion momenta given by 
2.8 <p n < 4.4 GeV In this range, a t0 + displays a soft 



decrease, which reflects itself in a soft increase of the nu- 
clear transparency. After including the effect of SRC (cal- 
culations referred to as RMSGA+SRC), the transparen- 
cies are about 2% larger for 12 C, 27 A1 and 63 Cu, and about 
4% larger for 197 Au. This reflects the fact that the medium 
effectively reduces the free pion-nucleon cross sections. 
It is obvious that the effect of SRC does not depend on 
the hard-scale parameter Q 2 . The CT mechanism, on the 
other hand, shows a strong Q 2 dependence with CT-related 
enhancements up to 20% at the highest energies. These 
calculations including CT are in very good agreement with 
the experimental data. The results overestimate the Au 
data somewhat, but the slope is in agreement with the 



0.85 




Q- CeV 



Fig. 2: [Colour online] The Q 2 dependence of the nuclear transparency for the A[e,e'n + ) process in 12 C, 27 A1, 63 Cu and 
197 Au. The dot-dashed curves are the results of the semi-classical model by Larson, Miller and Strikman 1 1 1511 . JLab data 
are taken from Ref. II16II . 



data. Upon comparing the RMSGA+CT results to the semi- 
classical calculations, we see that the slopes of both cal- 
culations are in excellent agreement, reflecting the use of 
the same quantum diffusion parametrization for the CT 
effect. There are some differences in the predicted values 
of the transparencies between the semiclassical and the 
RMSGA model. The RMSGA predictions are somewhat 
larger for the 12 C target, and evolve to smaller for the 
197 Au target. A more recent model developed by Kaskulov 
et al. II17II also finds excellent agreement between the cal- 
culations including CT and the data. 

Density dependence We can also apply the RMSGA model 
to map which density regions of the target nucleus are 
effectively probed in a certain nucleon-knockout reaction 
II18II . The scattering parameters entering in the Glauber 
profile function of Eq. ^ show little energy dependence 
for nucleon or pion momenta above 1 GeV This is re- 
flected, for instance, in the soft energy dependence of the 
regular RMSGA transparencies in Fig. [2J The stronger Q 2 
dependence observed in the data can be explained by the 
CT phenomenon. Together with the ability to deploy the 
RMSGA framework in a variety of reactions, this allows 
us to make statements about the role of nuclear attenu- 
ations on the effectively probed densities for a broad en- 
ergy range. 

For single-nucleon knockout reactions, we can study 
the density dependence of the reaction in a factorized ap- 
proach by calculating 5(r, 6>), which represents the con- 



tribution to the distorted momentum distribution around 
r and 6 and is defined as follows: 
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In this equation, <j) nKm (x) is the single-particle wave func- 
tion of the struck nucleon (with quantum numbers [riKni)), 
the missing momentum p m = q — p is defined as the dif- 
ference between the momentum transfer and the final nu- 
cleon momentum, and the function D(r,6) is defined as 



e -ipmr 
D{r,6)= I d<t>r 2 sm6— _u(p m , s )^ T (r)</v. m (r) . (9) 



(2ny 

In the factorized approach, the cross section is propor- 
tional to the distorted momentum distribution p° K (p m ), 
and <5(r, 0) thus provides a measure for the contributions 
to the cross section of the different density regions in the 
nucleus. 

Results for 5{r,Q) for the 12 C(e,e'p) and 12 C(p,2p) 
reactions are shown in Fig. [3] with 9 measured from the 
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Fig. 3: [Colour online] The function 5(r,6) for the 12 C(e,e'p) and 12 C(p,2p) reactions. The energy transfer for both 
reactions is 1.5 GeV and the three-momentum transfer q is tuned to probe the maximum of the momentum distribution 
(i.e. p m =0 MeV for knockout from the sl/2-orbit and p m =115 MeV for removal from the p3/2-orbit). The proton is 
ejected along q for the (e,e'p) results. For the (p,2p) calculations, both protons in the final state have 1.5 GeV kinetic 
energy and are detected under an angle of 32° on opposite sides of the incoming proton. This particular kinematic 
situation is often referred to as coplanar and symmetric kinematics. For the sake of reference, the proton root-mean- 



square radius in C as determined from elastic electron scattering is (r 



1/2 



:2.464±0.012fm[19]. 



direction of momentum transfer and r denoting the dis- 
tance relative to the center of the target nucleus. The 
energy transfer is 1.5 GeV and the polar angle of the 
ejected protons are determined such that the kinematics 
probes the maximum of the undisturbed momentum dis- 
tribution. In the relativistic plain-wave (RPWIA) limit, no 
ISI or FSI are present (^(r) = 1) and <5(r, 0) is equal for 
(e,e'p) and (p,2p). For the selected kinematics described 
in the caption to Fig. 3, there is a symmetry-axis for 6 = 
90° . In the absence of any effect stemming from nuclear 
attenuation the forward and backward hemispheres equally 
contribute. After including the FSI and/or ISI in the RMSGA 
approach, we observe that the nuclear attenuations re- 
duce the value of 5{r,Q), shift the maximum values to 
higher values of r, and also induce an asymmetry in the 
9 direction. The biggest contributions to the cross section 
stem from the forward hemisphere. The relative contri- 
bution from the interior regions of large target-nucleus 
density gets reduced due to its sensitivity to strong FSI 
mechanisms. All these effects are more pronounced in 
the (p,2p) reaction, where three (one incoming, two out- 
going) protons are subject to the attenuations and the 
biggest contributions are close to the nuclear surface. 

To formulate the equivalent of Eq. §8§ for two-nucleon 
knockout reactions such as A{y,pp~) in a factorized ap- 
proach, we have to assume that the proton pair resides 



in a relative S -state. This is a reasonable approximation 
as investigations of the 16 0(e,e'pp) reaction at the elec- 
tron accelerators in Mainz [20,21] and Amsterdam 11221 
l23l|24ll have clearly shown that pairs of protons are solely 
subject to short-range correlations when they reside in a 
relative S state under conditions corresponding with rel- 
atively small cm. momenta P (or, the initial protons are 
very close and moving back-to-back) . This assumption al- 
lows us to write a distorted momentum distribution for 
the nucleon-nucleon pair: 



Pn lKl ,n 2 K 2 ^J 2-1 
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= dRd65<iR,6) , (10) 

with P (p) the center of mass momentum (relative mo- 
mentum) of the outgoing nucleon-nucleon pair. The miss- 
ing momentum P m = P— q can be interpreted in the quasi- 
free approximation as the center of mass momentum of 
the correlated pair before the photon interaction. 

Fig.[4|shows 5(R,6) for the dual proton knockout on 
12 C, with proton kinetic energies of 1.5 GeV in the final 
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Fig. 4: [Colour online] The function 5(R,9) for the ex- 
clusive 12 C(y,pp) cross section. We consider an energy 
transfer of 3 GeV and a three-momentum transfer q that 
is tuned to probe the maximum of the momentum dis- 
tribution Pn 1 K 1 ,n 2 K 2 (P) (i- e - P = ® MeV for knockout from 
the (sl/2-sl/2)- and (p3/2-p3/2)-orbits). We consider 
coplanar and symmetric kinematics, i.e. the two escaping 
protons have the same energy and polar angle pq , but 
escape from the opposite side of q 



state. We can again quantity the effect of FSI by compar- 
ing the RMSGA to the RPWIA results. It is clear that the 
biggest values of 5(R,9) are situated much closer to the 
center than for the one-nucleon knockout reactions. Con- 
sequently, the A(y,pp) reaction succeeds in probing the 
high density regions of the target nucleus. 



4 Conclusion 

We have shown a selection of results obtained in a model 
based on relativistic multiple-scattering Glauber scatter- 
ing theory. The model has no free parameters and can be 
applied to variety of reactions, with leptonic and hadronic 
beams and outgoing nucleons and/or pions. Our calcula- 
tions show that relativity plays a rather modest role in 
the magnitude of nuclear attenuation. We showed very 
good agreement of pion transparencies with calculations 
including colour transparency. Secondly, we exploited the 
robustness of the model to explore which target-nucleus 
densities can be effectively probed in knockout reactions 
involving one, two and three protons. We find that the 
A{j,pp) reaction probes the interior of the target nucleus, 
the A(p,2p) is rather peripheral, whereas the A(e,e'p) is 
somewhat intermediate between these two. 
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